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Abstract 

The inverse scattering method is applied to the investigation of 
the equilibrium configuration of black holes. A study of the boundary 
problem corresponding to this configuration shows that any axially 
symmetric, stationary solution of the Einstein equations with discon- 
nected event horizon must belong to the class of Belinskii-Zakharov 
solutions. Relationships between the angular momenta and angular 
velocities of black holes are derived. 



1 Introduction 

In the present paper, we study solutions to the Einstein equations in a vac- 
uum that are stationary, axially symmetric, and asymptotically flat. The 
main purpose is to elaborate a description in which two classes of ideas, black 
hole theory and the theory of completely integrable equations, are unified. 
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In black hole physics, the well-known result claims that "black holes have 
no hair". In other words, each black hole is uniquely determined by its 
mass and angular momentum and is nothing but the Kerr solution. Many 
authors contributed to this result and we do not have room to give the 
exhaustive references. For our purpose, only the final stage of the proof is 
important, namely, Carter's classification of all axially symmetric solutions 
having connected horizon fl]]. Carter demonstrated that these solutions must 
solve the boundary problem for a system of elliptic nonlinear equations, the 
Kerr solution being one of the possible solution to this boundary problem. 
The uniqueness of this solution was proven in ||. 

Can one do without requiring the connectedness of the horizon? In other 
words, do solutions exists that correspond to an equilibrium configuration of 
black holes? To the best of the author knowledge, these questions still have 
no exhaustive answers (though some results can be found in 0, M). Also, 
these questions can be formulated in terms of the boundary problem from 
[H] because, there, the connectedness of the horizon played no role. This 
boundary problem is formulated in the following section. 

It is well known that the Einstein equations with two commuting symme- 
tries belong to a wide class of systems that can be integrated by the methods 
of the inverse scattering problem. This was shown in || by Belinskii and 
Zakharov, who also investigated the axially symmetric case for which the 
2iV-soliton solution was constructed |J. This solution was interpreted as 
the solution corresponding to the N Kerr-NUT black holes. As we show 
in the present paper, the 2iV-soliton solution of Belinskii-Zakharov indeed 
contains all possible solutions (if any exist) corresponding to an equilibrium 
configuration of rotating black holes. Therefore, we reduce the question of 
the existence of solutions with disconnected horizon to the investigation of 
some subclass of the soliton solutions. This subclass is parameterized by 
the distances between black holes, the angular momenta, and the masses of 
the black holes. However, in general case, solutions from this have a conical 
singularity on the symmetry axis, which hinders the existence of solutions 
with disconnected horizon. It is still unclear whether it possible to choose 
the parameters in a way that removes this conical singularity. Most likely, 
the answer is negative @, |J. 

The presence of the conical singularity does not make these solutions 
physically meaningless. On the contrary, the conical singularity itself has 
the physical sense of the interaction force between the black holes. This 
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interpretation was proposed by Weyl; for details see |7]]. 

2 Boundary conditions 

In this section, we present some basic facts about the stationary axially 
symmetric solutions to the Einstein equations and formulate the boundary 
problem corresponding to an equilibrium configuration for black holes. The 
results of |1J are crucial for us and we refer the reader there for details. 

Recall that Lorentz manifold is stationary and axially symmetric if it 
possesses two commuting one-parameter isometry groups that are isomorphic 
to R and SO (2) respectively. In other words, there exist two Killing vector 
fields k a and m a , which commute with each other, such that the vector k a 
is time-like and the vector m a is space-like. For the stationary and axially 
symmetric vacuum (or electrovacuum) space-time one can always choose the 
coordinate system in which the metric acquires the following form: 

ds 2 = -Vdt 2 + 2Wdtd<p + Xd<p 2 + ^(dp 2 + dz 2 ). 

P 2 

Here the metric coefficients depend only on p and z, and k a d/dx a = d/dt, 
m a d/dx a = d/dcj). Then V, W,X,p have a geometric meaning, 

V = -k a k a , X = m a m a} W = k a m a , (2.1a) 

VX + W 2 = p 2 = -p ab p ab (p ab = 2k [a m b] ). (2.16) 

The multitudes of the Killing vectors are normalized as follows: V — > 1 at 
infinity and X \x ' a ~ * ^ on ^ ne symmetry axis. The coordinates 0, p and z 
become cylindrical (Weyl) coordinates. 

The set of points with p = 0, X = is the symmetry axis, while the set 
of points with p = 0, X > is the event horizon. Let the event horizon have 
iV connected components and l a be an isotropic vector that is orthogonal to 
the event horizon. Then, in each connected component of the horizon, we 
can choose such a normalization of l a that 

i a = k a + n im a , 

where Qi is some constant whose physical meaning is the angular velocity of 
the black hole. Let z\,Z2, . . . , z^n be the points of intersection of the horizon 
and the symmetry axis enumerated in increasing order. 
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We pass to a new coordinate system, 

p 2 = (\ 2 -m 2 )(l-p 2 ),m. 

Z2i + Z2i-1 . 

z = \p. 

In this coordinate system, the necessary conditions of regularity of the sym- 
metry axis and the horizon are formulated as follows: 

X(\,p) = (l-p 2 )X(\,p) 

W^X,p) = (X 2 -m 2 )(l-p 2 )W(X,p) (2.2) 
V^\,p) = (\ 2 -m 2 )V(\,p) 
Here X, W, V are smooth functions nowhere equal to zero and 

= -l a l a , = l a m a . 

From (2.2) one can easily obtain the boundary conditions for the first group 
of Einstein equations, 

d*pdgg- 1 = 0, g = ( ^ ^ , (2.3) 

where * is the Hodge operator *dp = dz, *dz = —dp. Indeed using (2.2), one 
can easily prove that 

1 f0 0(1) 



P9, P 9 = { Q g ) ' ZE r ' (2 ' 4a) 

Qipg^g- 1 ^ 1 -(oil) ) ' p ^ z e 4 ^ = ( 1*) ' (2 ' 46) 
pg^g- 1 = 0{l), p^0,^£ (2.4c) 
Here T is the symmetry axis consisting of iV + 1 connected components, 

T = \jT J = R\\Jl t , j = l,...,N + l, i = l,...,N, I i = (z 2i - 1 ,z 2i ). 

The symbol 0(1) denotes a uniformly bounded function on corresponding 
interval. It follows from (2.2) that (2.4c) tends to zero almost everywhere 
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except the points Zk', however for our purposes, a uniform boundness suffices. 
The function g(z,p) is taken to be smooth at all points except the points 
(*fc,0). 

As we demonstrated below, Eqs. (2.4) completely determine the solution 
to Eqs. (2.3). At the same time fli and z$ are the independent parameters 
of the boundary problem; these parameters can be choose arbitrary. The 
behavior of g at infinity is discussed at the end of this section. 

An alternative approach, in which the main parameters are angular mo- 
menta rather than angular velocities, exists 0, ||. To show this let us intro- 
duce the Ernst potentials 

P * dgg- 1 = ( ~ dY X ) , d(Y 21 - Y 12 ) = 2dz. (2.5) 



-dY dY 2 . 

Here Y is the Ernst potential that is determined by the space-like Killing 
vector field, while Y is the Ernst potential determined by the time-like Killing 
vector field. Then system (2.3) can be rewritten in equivalent form, 

<i(^)-£*aTAaT = 0, 



d('^-) + ^*dXAdY = 0, (2.6a) 
JL X 



p * dY 
~X 

From (2.4a) or (2.2) we can see that 



dQ = - v2 , W = QX. (2.66) 



pd p \nX^2, p^O, zeT; Y\ Fi = q, (2.7) 

where q are some constants that are independent parameters. In || [7|, |], 
it was proved that (2.6) has a unique solution satisfying (2.7) and some 
condition at infinity for all Z{ and q. In fact, the latter condition is equivalent 
to the asymptotic flatness of the metric. Note that (2.4b) follows from (2.6b) 
provided that fit is defined as follows: 

fik = -a, (2.8) 
Thus, Qj is a function of zt and c,-. 
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The constants q unambiguously determine the angular momentum of all 
black holes. Indeed, let us define the angular momenta of a black hole using 
the Komar form: 

Li = — p / m a *dS ah} 

l07T JSi 

where Si is a two-surface surrounding the black hole. Choosing Si to be the 
surface of revolution of the curve connecting the components of the axis 
T i+ i and Tj, we obtain 

L t = lf -*(XdW- WdX) = \ ! dY = l(c i+1 -a). (2.9) 

o Jd p o Jd o 

The mass of the black hole can be defined as the following integral I]: 

Mi = -±- [ k a ' b dS ab . 
Sir J Si 

Proceeding as for finding the angular momentum, we obtain that 

Mi = j[ - * {XdV + WdW) = ~f dY 12 . 

Contracting the contour C, to the horizon and using (2.4b), we find that 

Mi = j f (2 + QiY z )dz = mi + 2QiL u (2.10) 

where mi = (z 2j - z 2j -i)/2. 

At infinity, we impose the following conditions: 

W = p 2 0(- 3 ), X = p\l +0(-)). (2.11) 

where r = y/p 2 + z 2 . Asymptotic formulas (2.11) are assumed to be differ- 
entiable at least twice. Formulas (2.11) mean that the metric tensor g tends 
to the Minkowski tensor in cylindrical coordinates. From (2.11) we obtain 




(2.12a) 
(2.126) 
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When determining the angular velocities from (2.8), we should normalize Q 
in accordance with (2.11), i.e. Q = 0(l/r 3 ). 

The second group of Einstein equations allows one to determine the coef- 
ficient e@ from the matrix g. Using (2.4) or (2.7), one can show that d z (3 = 
forp = and zGT, i.e. /3|r ; = &i, where 6j are some constant. The conical 
singularity on the symmetry axis is absent iff 6^ = 0. However, hi cannot 
be treated as independent parameters since they are functions of z^ and c^. 
In the preset paper, we restrict ourselves to the study of boundary problem 
(2.4), (2.12) and do not discuss the properties of &j. 

3 Auxiliary linear problem 

System of equations (2.3) is the compatibility condition for the following pair 
of matrix linear differential equations 0, |(J: 

Diip = — ip, D 2 ip = — — i). 3.1 

Here D\ and D 2 are the commuting differential operators: 

Di = d z ^ „ <9^, D 2 = d p + 'l UJp - d m 

and uj is a complex parameter that does not depend on the coordinates. We 
also use the U — V pair representation in which a; is a dependent parameter. 
To be more precise, let uj be a root of the equation 

uj 2 - 2uj(k - z) - p 2 = 0, (3.2) 

where k, in turn, is independent spectral parameter. Using (3.2), one can 
easily check that 

UJ 2 + p 2 U> 2 + p 2 

Passing from ipiuj) to ip'ik) = ip{uj{k)), we obtain from (3.1) that 

d z m = A(z,p,k)m, A= p2g ^\~ j UJp 2 g ^\ (3.4a) 
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d P m = B(z, P ,k)m, b = ^4t?^ ( 3 - 4& ) 

Hereafter, we omit the prime for brevity. It is worth mentioning that Eqs. 
(3.4) are equivalent to Eqs. (3.1) only if we take u> to be the multivalued 
function in (3.4). Fixing the branch of the root in (3.4), we find the solution 
to system (3.1) only in the analyticity domain of uj(k). 

In the present paper, we follow the general scheme for investigating inte- 
grable equations ||. 

Since Eq.(3.2) is invariant with respect to the transformation u — > —p 2 /u), 
we can fix the branch of the multivalued function u>(k), stipulating that the 
inequality \u\> p holds. Then, from (3.2), we obtain 

u — > 2{k — z), p — > 0; uj — > 2{k — z), z — > oo; u> — > 2{k — z), k — * oo (3.5) 

After choosing the branch of the root, we can introduce the monodromy 
matrix T(z, y), which, by definition, is a solution to (3.4a) such that T(y, y) = 
I. Note that for p — > 0, 

Here we took into account (3.5) and boundary condition (2.4). Hence Eq. 
(3.4a) can be easily integrated at p = 0. As a result, the explicit formulas 
for the monodromy matrix are 

/, Y(z)-Y(y) \ 

T(z,y)=( g-v) , z,yeT k , (3.7) 

V U k-y ' 

where T k is the connected component of the symmetry axis and 

( ¥r 0\ . 
T(z,y)=a t Y(z)-hy) z,yeli. (3.8) 

V 2(k-y) J 

Let e(z, p, k) be the solution to (3.4) with the Minkowski tensor in cylin- 
drical coordinates (V = 1,W = 0,X = p 2 ), 

e(z,k) = (l w( ° Jfc) ) i a,e = A,e,^o=(o __!_)• (3-9) 
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Let us define the Jost functions and reduced monodromy matrix, 

^{z,k) = lim T(z,y)e(y), (3,10) 

y— >±oo 

We do not reproduce the explicit dependence on p. The functions \E' ± satisfy 
the following integral equations: 

ty-(z) = e(z)+ I* e(z)e- 1 (x)A'(x)m-(x)dx, (3.12) 

J— CO 
/•CO 

^ + (z) = e(z) - J e(z)e- 1 (x)A'(x)^ + (x)dx, (3.13) 

where 

m =A { z ) -Mz) = ^--^ ? ( P9 ,^-§ »)). (3.14) 

Limit (3.11) exists at least for |Im/c| > p. Further, using (3.4b), one can 
show that T(k) does not depend on p. The basic property of the monodromy 
matrix reads 

T(z, y) = T(z, z 2N )T(z 2N , z 2N -i) ■ ■ ■ T(z u y), z E T N+1 , y e T 1 , (3.15) 

where T^+i and Ti are the extreme components of the symmetry axis. Then, 
from (3.15), (3.11), (3.8), (3.7) and (3.5) we get that 

T(k)=D N+1 [] TjDj. (3.16) 

j=l,...,N 

Here „ M . 

r-( - iM ' n '\ ft - n -Bi 

V (fc~Z2j)(fc — Z2j-l) k — Z2j ' 

where the constant Dj are defined as follows: 

D N+1 = (Y(oo) - Y(z 2N )), Dj = (y(^-i) - Y(z 2j _ 2 )), 
D 1 = (Y(z 1 )-Y(-oo)) 
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In (3.16), we took into account identities (2.9) and (2.10). Notice also that 
detT(ife) = 1. 

The cut of u(k,z,p) is the segment that connect the points z + ip and 
z — ip. Therefore, \1/ ± (A;) are analytic functions in k as \lmk\ > p and 

tf-(jfe) = V+(k)T(k). (3.17) 

The function ty + (k) (ty~(k)) can be analytically continued in the domains 
Rek < z (Rek > z). Using (3.5), we can see that for k — > oo, 

^ ± (A;)e- 1 (A;) -> / (3.18) 

Substituting k = z + (w 2 — p 2 )/2u for fc, we pass from \l/ ± (/c) to ^ ± (cu). 
Then the function ^/ ± (a;) become solutions to Eqs. (3.1) in the domain 
\u\ > p, while is analytic in uj as Rew < 0, \u\ > p and ^~{uj) is 

analytic as Reuo > 0,\uj\ > p. Furthermore, it follows from (3.18) that 

^ ± (cu)e~ 1 (cj) -> / (3.19) 

as w ^ oo. Though \l/ ± (u;) are determined for < p as well, they do not 
satisfy system (3.1) in this domain. Therefore, our next aim is to continue 
ty ± (uj) into the domain \uj\ < p in a manner that preserves Eqs. (3.1). 

Let iVi(k), 0J 2 {k) be the roots of Eq. (3.2) and let Reui(k) < and 
Reu 2 (k) > 0. Then the cuts of coi j2 (k) are half-lines going from points z + ip, 
z — ip to infinity in a direction that is perpendicular to the real axis. Hence, 
the functions ooi i2 (k) are analytic for |ImA;] < p. Note that oo\(k) = u(k) 
for ReA; < z and ou 2 (k) = c^(fc) for Rek > z, whence the functions \l/ ± (/c) 
are continued analytically into the strip \lmk\ < p. Further, let ^i(k) and 
^ 2 (k) are the solutions to Eqs. (3.12) (with u\ substituted for a;) and (3.13) 
(with u 2 substituted for u), respectively . 

For \Imk\ < p, the functions ^ + (£;) and ^i(k) (^~(k) and ^ 2 (k)) are 
the solutions of the same differential equation (3.4a). Hence, 

*i{k) = * + (k)T 1 (k), *-(k) = *t(k)T2(k), (3.20) 

where the matrices T 1>2 (k) do not depend on z. Since the solutions to the 
integral equations (3.12) and (3.13) automatically satisfy (3.4b) (this follows 
from boundary conditions (2.12) and the fact that e(z,p,k) is a common 
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solution of Eq.(3.4) with the Minkowski tensor), we conclude that T\^{k) do 
not depends on p as well. Moreover, at p — > oo, 

u>i(k) -> -p + (fc - z) + O(-), cu 2 (A;) -> p + (fc - z) + O(-). 

■ p p 

Therefore, accounting for Eq.(2.12), we derive that 

lim e^ l (k)^(k) = lim er^AOtf+ffc) = /, 

p— >oo p^oo 

lim e^(k)^!~{k) = lim \k)^t(k) = I, (3.21) 

and, hence, Ti(k) = T 2 (k) = I. Here e\{k) = e(u>i(k)) and C2(k) = e(u>2(k)). 
In other words, for |ImA;| < p, we obtain 

tf-(jfe) = tf-(jfe) = tf+(jfc). (3.22) 

The function ^(/c) and ^(A;) continued analytically into the domains 
Kek > z and Re/c < z respectively. Hence, the function ^^(u) is analytic for 
\uj\ < p, Keu < 0, and the function ^(u;) * s analytic for \u\ < p, Reu; > 0. 
As these functions are the solutions to Eq.(3.1) in the domain Reu; < 
or in the domain Reu> > (^f), it follows from (3.22) that \l/ + (u;) can be 
analytically continued into the half-plane Reu; < and remains a solution of 
(3.1), and fy~(ou) can be analytically continued into the half-plane Reu; > 0, 
and also remains a solution of (3.1). 

System (3.1) is invariant w. r. t. the transformation 

*(z,p,u)^g*-\z,p,-£), (3.23) 

which is valid because the matrix g is symmetric (the tilde in (3.23) denotes 
transposition). Reduction (3.23) means that ^~(k) and g[ty±(k)}~ 1 are the 
solutions of the compatible pair of equations (3.4) with uj — u> 1 . However, at 
P -> oo, 

e^(k)ge^(k) - -/. (3.24) 

Then from (3.21) and (3.24) we have that ^~(k) = -^[^(A;)] -1 ,or, equiv- 
alent ly, 

¥» = (3-25) 

uj 
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The functions and \& (uj) satisfy the compatible system of equa- 

tions (3.1). Therefore, the combination [^/ + (uj)]^ 1 ^/^ (uj) depends only on 
k = z + (to 2 - p 2 )/2lu (Dik = D 2 k = 0). Then, by virtue of identity (3.17), 

= V + (uj)T(k). (3.26) 

Equations (3.26) and (3.25) show that 

T(k) = f(k). (3.27) 

The monodromy matrix, T(k) depends on 3N + 1 parameters except z^. 
We treat equality (3.27) as the system of 2iV+l nonlinear algebraic equations 
for the constants Dj, Lj and Qf 

JV+l N 

J2 Dj + J2^jMj = 0, Res 2fc T 12 (A;) = Res Zfe T 21 (A;). (3.28) 

3=1 3=1 

Assume Dj can be excluded from (3.28); then the remaining N equations 
give us the connection between the angular velocities and angular momenta. 
For instance, for the case of a single black hole, it follows from (3.28) that 

D 1 = D 2 = -2MiHi, fti = t— — -. (3.29) 

The expression for Vt\ in (3.29) is known; it establishes the connection be- 
tween the angular velocity and angular momentum of the Kerr black hole. 
The relationship is expressed by 

r^L ^ ^1 7i t2 2 2 

^1 = 777 76 q> a = TT-> M 1 = m 1 + a . 

(Mi+mi) 2 + a 2 Mi 1 1 

As in (3.29), equality (2.10) is also taken into account here. Then the mon- 
odromy matrix reads 



T(k) 



I _ 2M i _|_ 2Mi(Mi-mi) 2Li 

k— z\ (k-z 1 ){k-z 2 ) (k—zi)(k-z 2 ) 

2Li y + 2M i _|_ 2Mi(Mi-mi) 

(fe— zi)(fc— Z2) k-z 2 (fe— Zl)(fe— Z2) 



Let us summarize the results of this section. Let boundary problem (2.4), 
(2.12) have a solution. Then there exists a piecewize analytic matrix ^(u) 
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($((*)) = ^ + (uj), Reu < and ty(u) = *~(o;), Reus > 0) that satisfies the 
compatible system of linear equations (3.1), the conjugation condition on the 
imaginary axis, 

\tj_(uj) = # + ( w )T(Jfe), Retu = (3.30) 
and the normalization condition at infinity, 




(3.31) 



In (3.30) #±(u;) = lim e ^ *(w T e) (e > 0). 

The only singularities of the matrix T(/c) as a function of are simple 
poles at the points tof = (zi — z) ± ^/ (z^ — -z) 2 + p 2 . Since the matrix 

(I X t JL) 

is regular at lg> = and tends to / as — > oo (which follows from explicit 
form of T(A;) and from Eqs. (3.27), (3.28)), we conclude that 

where do not depend on u;. 

Formula (3.32) demonstrate that a rational-in-cu solution to the auxiliary 
linear problem corresponds to each solution with a disconnected horizon and, 
hence, each such solution should be contained in the Belinskii-Zakharov class 
of solutions ||. 

Assume that for any Qi,Zi (or Li,Zi) the system of nonlinear equations 
(3.28) has a unique solution. Then a unique symmetric matrix T(k) cor- 
responds to each solution of the boundary problem (2.3), (2.4), and (2.12). 
However, since the solution to the Riemann problem (3.30), (3.31) is unique 
and g is unambiguously reconstructed by ^(co) (see 3.25), we conclude that if 
a solution to the problem (2.3), (2.4), (2.12) exists, then it is unique as well. 
In particular, for the case where only one black hole is present, we obtain a 
new proof of the uniqueness of the Kerr solution. 
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